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Abstract 

The problem of errors, arising due to finite BPM resolution, in the difference 
orbit parameters, which are found as a least squares fit to the BPM data, is one 
of the standard and important problems of accelerator physics. Even so for the 
case of transversely uncoupled motion the covariance matrix of reconstruction 
errors can be calculated "by hand", the direct usage of obtained solution, as 
a tool for designing of a "good measurement system", does not look to be 
fairly straightforward. It seems that a better understanding of the nature of 
the problem is still desirable. We make a step in this direction introducing 
dynamic into this problem, which at the first glance seems to be static. We 
consider a virtual beam consisting of virtual particles obtained as a result 
of application of reconstruction procedure to "all possible values" of BPM 
reading errors. This beam propagates along the beam line according to the 
same rules as any real beam and has all beam dynamical characteristics, such 
as emittances, energy spread, dispersions, betatron functions and etc. All 
these values become the properties of the BPM measurement system. One can 
compare two BPM systems comparing their error emittances and rms error 
energy spreads, or, for a given measurement system, one can achieve needed 
balance between coordinate and momentum reconstruction errors by matching 
the error betatron functions in the point of interest to the desired values. 
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1 Introduction 



The determination of variations in the transverse beam position and in the beam 
energy using readings of beam position monitors (BPMs) is one of the standard 
and important problems of accelerator physics. If the optical model of the beam 
line and BPM resolutions are known, the typical choice is to let jitter parameters 
be a solution of the weighted linear least squares problem. Even so for the case of 
transversely uncoupled motion this least squares problem can be solved "by hand", 
the direct usage of obtained analytical solution as a tool for designing of a "good 
measurement system" does not look to be fairly straightforward. It seems that a 
better understanding of the nature of the problem is still desirable. 

A step in this direction was made in the paper [1] , where dynamic was introduced 
into this problem which in the beginning seemed to be static. When one changes the 
position of the reconstruction point, the estimate of the jitter parameters propagates 
along the beam line exactly as a particle trajectory and it becomes possible (for 
every fixed jitter values) to consider a virtual beam consisting from virtual particles 
obtained as a result of application of least squares reconstruction procedure to "all 
possible values" of BPM reading errors. The dynamics of the centroid of this beam 
coincides with the dynamics of the true difference orbit and the covariance matrix 
of the jitter reconstruction errors can be treated as the matrix of the second central 
moments of this virtual beam distribution. 

In accelerator physics a beam is characterized by its emittances, energy spread, 
dispersions, betatron functions and etc. All these values immediately become the 
properties of our BPM measurement system. From now one can compare two BPM 
systems comparing their error emittances and error energy spreads, or, for a given 
measurement system, one can achieve needed balance between coordinate and mo- 
mentum reconstruction errors by matching the error betatron functions in the point 
of interest to the desired values. 

This dynamical point of view on the BPM measurement system was explored 
in pp in application to the case of transversely uncoupled nondispersive beam mo- 
tion and in this paper we continue this study adding energy degree of freedom^ 
The paper by itself is organized as follows. In section 2 we introduce all needed 
notations, formulate the problem and give its standard least squares solution. As 
a new element, we formulate the necessary and sufficient conditions for the BPM 
system to be able to distinguish between transverse and energy jitters in terms of 
its three BPM subsystems. In section 3 (the core section of this paper) we make 
parametrization of the covariance matrix of the jitter reconstruction errors using the 
usual accelerator physics concepts of emittance, energy spread, dispersion and be- 

^It is clear, that such considerations, if needed, can also be done for the case of the fully coupled 
six dimensional motion. It is also clear that in similar fashion one can approach some other problems 
connected with the error propagation. It should not be necessary the BPM reading errors, it could 
be, for example, errors in the kick angles produced by the orbit feedback system. 
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tatron functions. We also show that the error dispersion is not simply one of the 
many dispersions which could propagate through our beam line. It, in analogy with 
the error betatron functions [1], is by itself solution of some minimization problem 
and is uniquely determined by transport matrices between BPM locations and by 
BPM resolutions. In section 4 we consider the measurement system which utilizes 
three beam position monitors (the minimum number of BPMs needed) and analyze 
in details effect of symmetries of the optics between BPM locations. In section 5 we 
continue the investigation of periodic measurement systems started in [1]. This time 
with the main accent on achievable energy resolution. And, finally, in section 6 we 
discuss application of the Courant- Snyder quadratic form as error estimator, even so 
in the case when energy degree of freedom is taken into account this quadratic form 
is not bound to be an invariant. 

2 Problem and Its Least Squares Solution 

Let us consider a magnetostatic beam line which is built from optical elements 
which are symmetric about the horizontal midplane y = 0. In such magnetic system 
the transverse particle motion is uncoupled in linear approximation, the vertical 
oscillations are dispersion free and errors in reconstruction of their parameters were 
already studied in [T], and in this paper we will examine together x-plane and energy 
degrees of freedom because they are connected through (linear) dispersion. 

We will use the variables z = (x, p, e)~^ for the description of the horizontal 
dispersive beam motion. Here, as usual, x is the horizontal particle coordinate, 
p is the horizontal canonical momentum scaled with the kinetic momentum of the 
reference particle and the variable e stays for the relative energy (or momentum) 
deviationo As orbit parameters we will understand values of x, p and e given 
in some predefined point in the beam line (reconstruction point with longitudinal 
position s = r) and as transverse and energy jitter in this point we will mean the 
difference 

Sz{r) = {6x{r), Sp{r), 6e{r))^ = (x(r) — x(r), p(r) — p{r), e — e)^ (1) 

between parameters of the instantaneous orbit and parameters of some predeter- 
mined reference (golden) trajectory (x, p, e)^. 

Let us assume that we have n BPMs in our beam line placed at positions 
si, . . . ,Sn and they deliver readings 



^The exact form of the variable e which we have in mind can be found in [2], but let us note 
that for the present study the particular form of this variable is unimportant. Let us also note that 
while in [1] the symbol e was used for the BPM reading errors, in this paper we prefer to use it for 
the relative energy deviation, and for the BPM reading errors we will introduce as new notation. 




(2) 
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for the current trajectory with previously recorded observations for the golden orbit 
being 



hi 



(3) 



Suppose that the difference between these readings can be represented in the form 

/ x{si) - x{si) \ 



Sb^ =^ be - bg 



\ 



X [ Sri. ) 3^ 



(4) 



where the random vector q = (q, . . . , has zero mean and positive definite 

covariance matrix Vc , i.e. that 



(?) = 0, v(?) = (?-?^>-(?>-(?>" 



K > 0. 



(5) 



The purpose of this paper is to study the influence of BPM reading errors ? on 
precision of reconstruction of jitter parameters under assumption that optical model 
of the beam line is known. The additional assumptions which we will make are: the 
covariance matrix stays constant and the BPM reading errors can be treated as 
independent from one measurement to the other. So BPM errors that are correlated 
from measurement to measurement (calibration and other systematic errors, drifting 
BPM readings and etc.) and fluctuations in BPM resolutions will be not considered. 
In practical applications these assumptions may or may not be realistic, but, first, 
they make the underlying mathematics almost trivialfl and, second, their satisfaction 
is, in some sense, one of the goals for the BPM and BPM electronics designers. 

Let Am{r) be a transfer matrix from location of the reconstruction point to the 
m-th BPM location 



A. 



amir) Cmir) gm{r) 

emir) dm{r) fm{r) 
n n 1 



am{r) dm{r) - Cm{r) e„,(r) 



(6) 



and let us assume that the Cholesky factorization 
matrix is known. As usual, we will find an estimate 



RJ of the covariance 



6z^{r) 



{Sx,{r), Sp,{r), Se,{r)) 



(7) 



for the difference orbit parameters ([T]) in the presence of BPM reading errors by 
solving the following weighted linear least squares problem 



■^Under these assumptions errors in the reconstruction process can be modeled as a sequence 
of independent identically distributed random variables (like in coin tossing) and therefore all 
probabilistic characteristics can be obtained studying errors in reconstruction of the result of only 
one measurement, but for all possible values of 
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mm 



M,-5z, - -Sb, 



(8) 



Here 



denotes the Euclidean vector norm, = R '^M and 



M 



/ ai(r) ci(r) gi{r) \ 
\ anir) c„(r) gn{r) J 



(9) 



The problem ([8]) always has at least one solution and, if we will assume that the 
matrix has full column rank (rank(Mj) = 3), then the solution of this problem 
is unique and is given by the well known formula 



6z,{r) = {Mj{r)M,{r)) ^ MJ (r) R;^ ■ 6 b,, 



= {M^{r)V-^M{r)) M"^(r) V"^ ■ 5 6, . (10) 

The calculation of the covariance matrix of the errors of this estimate (object of our 
main interest) is also standard and gives the following result 

V,{r) = ViSz,{r)) = {mJ {r)M,{r)y' = {M^ (r) V~' M (r))''. (11) 

Let us discuss in more details the important condition for the matrix to have 
full column rank. This condition will allow us to separate betatron and dispersion 
oscillations at the BPM locations and, therefore, will make our system applicable for 
measuring transverse and energy jitter. 

Because the matrix R, is nondegenerated, the rank of the matrix is always 
equal to the rank of the matrix M, and the matrix M, in the next turn, will have full 
column rank if and only if the Gram determinant F (a, c, g) of its column vectors 

a= (ai, . . . ,a„)^ , c= (ci, . . . , c„)^ , g = {gi,...,gn)^ (12) 

is not equal to zero. 

To find desired expression for the Gram determinant let us introduce Bjnk - 
transport matrix from the location of the BPM with index m to the location of the 
BPM with index k 

Bmk = ^fe 

With these notations and using Binet-Cauchy formula one can obtain after some 
straightforward manipulations 



<*ii 

*21 



*22 




(13) 
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r (a, c, g) = det (M^M) 



27 7^ 

*52 *12 



ij jk 
*12 *16 



l<i<j<fc<n 



1 " 2 1 

— i \^ /^a3*i ff.i'' _ ^.J'^^ — _ \^ IfJi _ (^A\ 

~ g 2-^ \^*52 *12 *12 *16y ~ g 2-^ l,*12 *16 *16 <*12y ) l-l-^J 

i,j,k=l i,j,k=l 

where 8852 (in the framework of the usual 6 by 6 matrix formahsm for the hnear 
beam dynamics) is the coefficient that connects variation of the particle path length 
with variation of the particle transverse momentum and which can be expressed 
using elements of the matrix Bij as follows 

*i ij ij ij /I cN 

<*^52 — *22 *16 <*^12 *26 ■ K^'^J 

From ( lT4l) one sees, that the matrix M will have the full column rank if and only 
if there are at least three beam position monitors with indices i, j and k such that 
the transport matrices between them satisfy the condition 

^52^12 ~ ^12^16 7^ (16) 

or (equivalently) the condition 

^12 ^16 ~ ^16 ^12 7^ 0- (1"^) 

Note that both conditions, f ll6p and f ll7p . involve elements of two transfer ma- 
trices, but while f lT6|) uses matrices between neighboring BPMs {Bij and Bjk), con- 
dition f lT7|) operates with the transport matrices from first to two remaining BPMs 
[Bij and -Bjfc). In simple words the condition f|T7|) . for example, means that one 
can not vary particle transverse momentum and particle energy at the first BPM 
location in such a fashion that these variations are invisible at the two downstream 
BPMs. 



3 Beam Dynamical Parametrization of 

Covariance Matrix of Reconstruction Errors 

Let A{ri, r2) be a matrix which transport particle coordinates from the point 
with the longitudinal position s = ri to the point with the longitudinal position 

s = r2 

/ mil rriu ttliq \ 
-4(ri, r2) = 77721 "^22 , ttih 77^22 - "^12 "^21 = 1- (18) 
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Similar to [T], one can easily show that for any given value of ? the estimate of the 
difference orbit parameters propagates along the beam line exactly as particle 
trajectory 

6z,{r2) = A{ri, r2)-6z,{ri), (19) 

as one changes the position of the reconstruction point. So again we can consider 
a virtual beam consisting from virtual particles obtained as a result of application 
of formula (|TOl) to "all possible values" of the error vector ?. The dynamics of the 
centroid of this beam Szq coincides with the dynamics of the true difference orbit 

6zo{r) = {Sz,{r)) = 6z{r) , (20) 

and the error covariance matrix flTT]) can be treated as the matrix of the second 
central moments of this virtual beam distribution and satisfies the usual transport 
equation 

V;(r2) = A{n, ra) V,{n) A^(ri, ra). (21) 

Consequently, for the description of the propagation of the reconstruction errors 
along the beam line, one can use the accelerator physics notations and represent the 
error covariance matrix in the familiar form 



V, = {MjM,) ' 
















II 






^ "r 1 







-e, a, + A2 r/,,, r/^,, e, 7, + A^ r]^^ A^ rj^^^ | . (22) 
A'^.,. A^r/p,, A2 



As usual for the particle dynamics, this parametrization has two invariants (quan- 
tities which are independent from the position of the reconstruction point), namely 
transverse error emittance and rms error energy spread A^, which can be calcu- 
lated according to the formulas 



' A, = ./ , (23) 



where we have used the notations 

a, = R;^a, c, = Rfc, g, = Rfg (24) 
and r ("Ui, . . . Urn) is the Gram determinant of the vectors mi, . . . Um- 
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The error Twiss parameters, of course, remain the same as they were earher 
pubhshed in [1], namely 

P^r) = e, \\c,{r)\\l , a,{r) = e, (a,(r) ■ c,{r)) , 7,(r) = e, \\a,{r)\\l , (25) 
and for the new objects, the coordinate and momentum error dispersions, we have 
VxAr) = e,(^a,{r) ic,{r) ■ g,{r)) - /3,{r) (a,(r) ■ ^,(r)) ) , (26) 

VpAr) = ^,{oi,{r) {d,{r) ■ g,{r)) - 7,(r) (c;(r) ■ ^,(r)) (27) 

As it was shown in [1], the error Twiss parameters ( 125|) are not simply one of 
many betatron functions which could propagate through our beam line, they are by 
themselves solutions of some minimization problem and are uniquely determined by 
transport matrices between BPM locations and by BPM resolutions. And we would 
like to show, that the same is true also for the error dispersions fl26|) and fl27|) . 

Let rjxir) and ?7p(r) be some dispersions specified in the reconstruction point. 
Then the corresponding coordinate dispersion at the m-th BPM location can be 
calculated as follows 

Vx{sm) = am{r)ri^{r) + Cm{r) r]p{r) + gm{r). (28) 

Consider a vector 

D (r, ri^{r), r]p{r)) = R-~^ (?7x(si), . . . ,?7x(sn))'^ = Vx{r) + Vpir) + g^ (29) 
and a minimization problem 

niin D (r, r/^(r), r]p{r)) (30) 

By standard means it is not difficult to show that if F (a, c) 7^ then the solution 
of this minimization problem is unique and is given by the formulas fl2B]l and flTTI) . 

If, additionally, P (a, c, g) ^ then the minimum in (!30!l is bigger than zero 
(and is equal to zero otherwise) and the following identity holds 

2 1 

\\D{r, 7]^^,{r),7]p^,{r)) II2 = (31) 

Note that geometrically the vector rix{r) + rip[r) is nothing else as taken 
with an opposite sign projection of the vector g^ onto a linear subspace formed by 
vectors and and hence the vector D (r, r]x,<;{r), ripA^)) is orthogonal to both, 
vector and vector c^. 

To finish this section let us, for the case when readings of different BPMs are 
uncorrelated, i.e. when the covariance matrix is a positive diagonal matrix 
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K = diag(al^ al...,al) > 0, (32) 
write down the following useful expressions for the Gram determinants 



r(a„C„^J = - y ^523^12 35123516 =i ^ ^12^1^_^16^ ^(34) 

6,1^=1 V ^^^^-^^ / ^,riiiV ^^^^-^^ / 

which enter formulas fl23l) for the transverse error emittance and for the rms error 
energy spread. 

4 Three BPMs in Symmetric Arrangement 

Let us assume that we have three beam position monitors in our beam line which 
deliver uncorrelated readings with rms resolutions ai, a2 and (J3, and let B12 and 
-B23 be transfer matrices between first and second, and between second and third 
BPM locations respectively 

?^ii r-16 \ / mil "I12 mig 

B12 = I rai r22 , ^23= rragi | • (35) 

1/ \ 1 

When the phase advance between the first and the second BPMs or the phase 
advance between the second and the third BPMs is not multiple of 180° , i.e. when 

^?2 + ^ 0, (36) 

this system can be used for the measurement of the transverse orbit jitter with the 
transverse error emittance given by the following expression 

^^"^"^ . (37) 



(rfml^ + (y\ {■muri2 + r22mi2) + cr|r 



2 ^2 
12 



In order to be able to resolve both, transverse and energy, jitters simultaneously 
we have to assume, additionally to fl36l) . that 

mi2r52 - ri2mi6 7^ 0, (38) 

where the r^2 and r^i coefficients can be expressed using elements of the matrix 
B12 as follows 
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With fl36|) and fl38|l satisfied, we obtain for tlie square of the rms error energy spread 

^2 ^ ^1^12 + ^1 (^11^2 + ^22^12)^ + ajrl^ ^^^^ 
("^12^52 - runiief 

To complete description of the covariance matrix of the reconstruction errors 
fl22l) for the three BPM case, we also need formulas for the error coordinate and 
momentum dispersions, and for the error betatron functions. And although it is 
not very difficult to provide some formulas using (l25l) . (l26l) and ( 1271) . the results are 
not very informative and it is not easy to derive some nontrivial conclusions from 
them. So in this section we will give more digestible expressions for error dispersions 
and error betatron functions making an additional simplifying assumption about our 
measurement system that the transfer matrix B23 between the second and the third 
BPM is not an arbitrary beam transport matrix, but is obtained as a result of some 
symmetry manipulation with the transfer matrix between the first and the second 
BPM. 



4.1 Mirror Symmetric Optical System 

Let a magnet system between the second and the third BPMs be a mirror symmet- 
ric image of the magnet structure between the first and the second BPM locations. 
Then 



B23 = r2i rii -rsi . (41) 





The transverse error emittance of this measurement system is given by 

e = J ^1^2 ^3 .^2) 

|ri2| ^af + Aalrl^ + al ' 

and the error betatron functions at the BPM locations can be calculated as follows 

/X ■ ( ^ ^1 2a|r22 (1 + 2ri2r2i) + a^rn 

a,(si) = sign(ri2)--- ======= , 44 

^^2^3 + Aalrl^ + al 



= l^-l ■ ^ ■ / 2 1^1 ^ 2> (45) 
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a,{s2) = sign (ri2) ■ ( ^ - ^ ) ■ , , , 2 ' ^^^^ 



ftfe) = M--^- (47) 

/ N ■ f \ ^3 f^?rii + 2a|r22 (1 + 2ri2r2i) 
a^ss) = -sign (ri2) ■ ^======== . (48) 

If we will assume that BPM resolutions follow mirror symmetry of the system, 
which means that ai is equal to (T3, then, as it could be expected, the error Twiss 
parameters will satisfy the following symmetry relations 

/3?(s3) = /3?(si), a,{s3) = -af(si), ^^(52) = 0. (49) 
For the square of the error energy spread we have the following expression 



\ = 772 ' (50) 

^ '52 



and the coordinate and momentum error dispersions at the BPM locations are given 
below 

f s. _ (ri6 + 2ri2r5i) - 4 (t| ri2 r22 - 0-3^6 , . 

ri2- [O^l + 4(T2r22 + 
4 0"2 ^22 r^2 

"-'''■''^ = al + 4alrl, + al' ^^^^ 

. X (^friQ + 4 (T^ ri2 r22 r26 - erf (^e + 2ri2r^i) , . 

^P-^(^3) = .2 , A^2^2 , ^2x • (56) 
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Figure 1: Schematic layout of four bend chicane. 



One sees that if BPM resolutions will follow mirror symmetry of the system, they, 
similar to the error betatron functions, will satisfy 



(57) 



independently if is equal to zero or notfl One also sees that if mirror symmetric 
system can be used for energy jitter measurement (i.e. if 7^ 0), then the error 
dispersion is nonzero at the system entrance and exit, again independently if is 
equal to zero or not. 

As a more specific example, let us consider three BPMs integrated into four bend 
magnetic chicane, as shown by red circles in figure 1. For this system 



B 



12 




(58) 



where 



ri2 



Li + U + 



2Lo 



+ 



COs(v9) COS'^((/9) 



7^ 



and 



r52 



cos{(p) 



2L2 tan(v9/2) 



cos{(p) 



tan(y?) ^ 0. 



(59) 



(60) 



Therefore this system always can be used for the transverse and energy jitter mea- 
surement, and, as a concrete case, let us consider the first bunch compressor of the 
FLASH facility P H], which is the four bend chicane of the discussed layout. The 



*Let us remind, that the condition r26 = apphed in the symmetry point is the necessary and 
sufficient condition for the total transfer matrix of the mirror symmetric system to be achromatic. 
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2 3 4 

Energy BPM resolution (rms, microns) 



5 
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Figure 2: Resolutions of orbit and energy BPMs (shaded area) which are needed in 
order to be able to resolve energy jitter 5 ■ 10~^ in the first bunch compressor of the 
FLASH facility. BPMs are positioned as shown by red circles in figure 1. 

typical deflection angle for this chicane is about 18° and the distances L2 and L3 are 
equal each other and are equal to 0.5 m (see, for example, [5]). Let us assume that the 
first and the third BPMs (orbit BPMs) have the same rms resolutions ai = = aorb 
and for the second BPM (energy BPM) let us introduce the notation (T2 = o"enr- Let 
Arfes will be energy jitter resolution desired for the system. With these numbers and 
notations, and using the usual three sigma criterion (3A^ < A^es) we obtain from 
f l50|) the following inequality 



which gives us limitations on the range of the BPM resolutions which can provide 
the required precision for the energy jitter measurement. Figure 2, for example, 
shows the area of acceptable BPM resolutions defined by the inequality ( l6T]) for 



To finish the chicane discussion, let us move the first and the third BPMs into 
positions shown as green circles in figure 1. For this case 



al, + 2aL < 0.02663- A2 



'des ' 



(61) 



Ade, 



5 ■ 10-5. 




(62) 
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r22 = — = cos(v9), ri2 = — + L2 + L4 cos(v?), (63) 
rii cos^(v9j 



= -^2 tan(y?/2) - L4 sin(v9), r26 = -sin(v9), (64) 



r52 = L2 tan((y9/2) H tan(v?), (65) 

COS[(p) 

and one sees that this BPM positioning still can be used for the jitter measurement, 
because ri2 7^ and r^2 7^ 0, but both, the transverse error emittance and the error 
energy spread become larger (for the same BPM resolutions) than for the original 
BPM layout. Nevertheless, it is a good example of a mirror symmetric system for 
which r52 7^ rig and the total transfer matrix is not achromatic. 

4.2 Mirror Antisymmetric Case 

If a magnet system between the second and the third BPMs is a mirror antisym- 
metric image of the first part of the system, then 

r22 ^52 

B23 = I r2i rii rsi | . (66) 
1 

The transverse error emittance and the error beta functions remain the same 
as for the mirror symmetric case, but the measurement of the energy jitter is not 
possible anymore, independent of the BPM resolutions following symmetry of the 
system or not. The coordinate error dispersion is always zero at the BPM locations 
with the momentum error dispersion taking the values 

VpA^i) = VpA^s) = - — , VpAs2) = - — , (67) 

1^12 ^12 

which are independent from BPM resolutions. Note that this impossibility of the 
energy jitter measurement does not depend on the value of rig which could be zero 
or not|j 

4.3 Two Periodic System 

Let us assume that our measurement system is periodic, by which we mean 
that i?23 = Bi2- We named it in the title as two periodic owing the fact that two 
equal transfer matrices are involved, but, more correctly, it should be treated as 



^The condition rig = is the necessary and sufficient condition for tlic total transfer matrix of 
the mirror antisymmetric system to be achromatic. 
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a three cell system because we consider three BPMs. Note that general periodic 
measurement systems constructed from n identical cells will be studied in details in 
the next section, but with additional simplifying assumptions that the cell transfer 
matrix allows periodic beam transport and that all BPMs have the same resolutions. 

The transverse error emittance of the two periodic system can be expressed in 
the form 



1 criC2C"3 



ri2| + al-til{B^2) + (T. 



2 ' 
3 



(68) 



where 



tr, (5i2) = rn + , (69) 
and calculation of the error betatron functions gives the following results 

I . • / N ^1 a-| ■ tr^ (5i2) ■ (rn ■ tr^ (B12) - 1) + o-grii 
a,[si) = sign(ri2) , , o . ..^ . ' i^^) 



cr.aa ^aj + ■ tr^ (5,2) + ^ ^ 



a,[s2) = sign (ri2) ■ — , , ^ s ' l^^) 



^1^3 ^af + ai ■ til (5i2) + al 



as a\ + al ■ tr^ {B 



^1^2 ^al + ai ■ tr^ (B,,) + ^ ^ 



. s • ^ A ^3 (Ti r22 + (T2 ■ tr^. (5i2) ■ (r22 ■ tr^. (^12) - 1) , . 

adss) = -sign ri2 9 . 2 / n 2 " 

cTia2 v^a^ + ai ■ til (^12) + a| 

Let us assume in the following that BPM resolutions follow symmetry of the 
system, which, in the periodic case, naturally mean that ai = a2 = a^. In this 
situation /3^{si) and /^^(ss) are always equal to each other and, it seems, it is the 
only symmetry which does not require additional assumptions about coefficients of 
the cell transport matrix B12. 
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The error betatron functions will be cell periodic (will stay unchanged after trans- 
port through the first half of the system), if and only if 



and, if fl76l) is satisfied, then 



KiBu) = 1, (76) 



cos(/ip) = ±^ (77) 

and, therefore, 

sin(3/ip) = sin(/ip) ■ (4cos^(yUp) — l) =0, (78) 

where /ip is the cell phase advance corresponding to the periodic betatron function. 

The error betatron functions will be two cell periodic (will stay unchanged after 
transport through the whole system), if and only if 

tr^ (5i2) = tr, (5i2) , (79) 

which, when compared with ( !76|l . gives equation 

tr,(5i2) = 0, (80) 

as the condition for the "true two cell periodicity" (two cell periodic, but not one 
cell periodic). This condition means that the transverse part of the total system 
matrix (two by two submatrix located in the left upper corner) is equal to the 
minus identity matrix for which arbitrary incoming beta and alpha functions will be 
transported without changes through the system, but the error betatron functions 
will also bring the sum of the beta function at the BPM locations to the minimal 
possible value. 

To finish the discussion about error betatron functions let us note, that if in the 
matrix B12 the first two diagonal coefficients are to equal each other (rn = r22), 
then 

«f(s3) = -tt?(si), af(s2) = 0, (81) 

and one may say that in this situation the error betatron function becomes mirror 
symmetric. 

For the error energy spread and the error coordinate and momentum dispersions 
we have in the case of the two periodic measurement system the following formulas 



(^16 - ^^52) 



^2 ^ u,^u,-..^,^uj^u, ^ ^g2) 
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q-j (ri6 + rii (ri6 - rga)) + aj ■ tr^ (^12) • ((tr^. (^12) + 1) - rga) + ajne 

ri2 ■ (aj + al ■ tr^ (^12) + ^1) 



, , ^2-tr.(^i2)-(ri6 - r52) , . 

^'''^ '^ - al + cri ■ i.l (B,,) + al ' ^^^^ 



o-fri6 + cr| ■ tr^ (5i2) ■ (tr^ {B12) ■ r^2 + ?^22 (?^16 - ^^52)) + crlr->2 



(86) 



^^'^^ '^ a? + ■ tr^. (i?i2) + ai ' ^ ' 

VpA^s) = 

_ (ylr^2 + q-j • tr^ (^12) • ((tr^ (^12) + 1) r^2 - ri&) + er| {r^2 + ^"22 (^52 - ^le)) , . 

n2 ■ K + ai ■ tr^ (i?,2) + a|) '^''^ 

and, if resolutions of all three BPMs will be equal, the error dispersion will satisfy 
the equality 

VxA-^i) = VxA^s)- (89) 

The condition for the error dispersion to be cell periodic is more restrictive than 
for the error betatron functions, namely 

tr,(5i2) = -1, (90) 

and the condition for the "true two cell periodicity" is 

tr,(5i2) = -2, (91) 

which does not lead to any noticeable symmetry of the error betatron functions and 
which means that the transverse part of the cell matrix B12 is equal to the sum of 
the minus identity matrix plus some nilpotent matrix ( A^^ = ). 

Note that under condition fl90|) we have for the periodic cell phase advance the 
following relations 

cos(yLip) = sin(3/ip/2) = sin(/ip/2) • (2cos(/ip) + 1) = 0. (92) 
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4.4 Cell Followed by Switched Cell 

If a magnet system between the second and the third BPMs repeats the magnet 
system between the first and the second BPMs but with switched directions of dipole 
magnets, then 

' ru ri2 -ri6 
B23 = I -r26 | • (93) 

1 

In analogy with transition from mirror symmetric to mirror antisymmetric case, the 
transverse error emittance and the error betatron functions remain the same as for 
the two periodic measurement system, but, in contrast to mirror antisymmetric case, 
this system still can be used for the energy jitter measurement if 

ri6 + r52 ^ , (94) 

which, in particular, forbids the magnet system between the first and the second 
BPMs to be mirror symmetric by itself. 
For this measurement system we obtain 

(ri6 + rss) 

/ X (^I ine + ^52) , . 

_ cfI {fiQ - rii (ri6 + rsa)) + (t| ■ tr^ (^12) ■ ((tr^ (^12) - 1) rig - rgs) + crfrig 



ri2 ■ {(yf + al ■ til (^12) + f^l) 



(97) 



(T2 ■ (^12) ■ (ri6 + r52) , . 

VxAs2) = 2 ^ 2 /D X , T' 98) 



(^friQ - al ■ tr^ (^12) ■ (tr^ (^12) ■ r52 - r22 (rig + r^2)) - <ylr-,2 



ri2 ■ {(^1 + crl ■ ill (Bu) + (t. 



(99) 
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^ o"?^52 + (^2 ■ fag; (Bu) " ((tr^ (^12) - 1) ^52 - Tie) + (r52 - (r52 + rie)) 

n2 ■ (a? + ■ tr^ (B,,) + ai) ^ 

and one sees that for equal BPM resolutions the property 

VxAsi) = VxA^s) (102) 

is still preserved. 

There is no reasons to expect that coordinate error dispersion and simultaneously 
momentum error dispersion could stay constant at all three BPM locations (analogy 
of cell periodicity for the two cell measurement system) and, as one can check, there 
is no solution for that. Nevertheless, both error dispersions still can stay unchanged 
after transport through the whole system, if 

tr, (Bu) = 1 or tr, (^12) = 2. (103) 



5 Periodic Measurement Systems 

Let us consider a measurement system constructed from n identical cells assuming 
that the cell transfer matrix allows periodic beam transport with phase advance per 
cell /ip being not a multiple of 180°. Additionally, we will assume that BPMs placed 
in our beam line deliver uncorrelated readings, all with the same rms resolution abpm- 

Let us first consider the case when we have one BPM per cell (identically posi- 
tioned in all cells) with the periodic betatron function and the periodic dispersion 
function at the BPM locations equal to /3p(si) and 77a;,p(si) 7^ respectively. 

In this situation the formulas for the error transverse emittance and the error 
betatron function remain the same as was already published in [1], and the square 
of the error energy spread is given by the following expression 

A,^ = ■ Qnifip) , (104) 



where the function 



^ 1^ sin(?2/ip) 



Qni^^p) = , . . . rf^^— -—^ (105) 

^ ^ 1 sm{nfip) ^/ 1 sm(n/ip/2) 



n sin(/ip) \n sin(/ip/2) 

is defined only for n > 3|^ Note that for n > 3 this function could be extended by 
continuity for all fip not multiple of 360° where it becomes unboundedlll 

^For n = 1,2 the denominator in the formula (|105p is equal to zero independent of the value of 
the cell phase advance fXp 

^The nonnegative denominator in the formula ()105p is equal to zero not only when fip is a 
multiple of 360°, but also when n is even and, simultaneously, fip an odd multiple of 180°. 
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The coordinate and momentum error dispersions rjx^i; and rjp^^ at the BPM lo- 
cations are given below 



VxASk) = ?7x,p(Sl) ■ (l - UJnifip) ■ COS { "'+\ flp)) , (106) 



VpAsk) = VpAsi) - ^x,p(si)-^^y^(sin(^^±^/ip) - 



«p(si)-cos(^^±^/ip)) , (107) 



1 sin(n/ip/2)\ / 1 sin(nyUp; 



and one sees that while the coordinate error dispersion rj^^^ always have mirror 
symmetry 

VxA^k) = rjxA^n+i-k) , k = l,...,n, (109) 

the momentum error dispersion will be mirror antisymmetric 

Vpsi^k) = -Vpsi^n+i-k) , k = 1, ...,n (110) 

only in the case when ap(si) = and ?7p,p(si) = 0. 

Note, that the mean value of the coordinate error dispersion and the mean value 
of its squares satisfy the following relations 

The function ^„(/ip) is never smaller than one and is equal to one (reaches its 
minimum) only in the points 

,360° . oano\ 1 ri,...,ri-l if n is odd . . 

fip = k (mod 360 ), fc = < , n i n , i i -r • (112) 

'^^ n ^ ' 1^ 1, . . . , I - 1, 1 + 1, . . . ,n - 1 if WIS even ^ ' 

in which error dispersion coincides with periodic dispersion and which seem to be 
good candidates to be selected for improving resolution of the energy jitter measure- 
ment (see figure 3), if we are free in choosing the cell phase advance while, for some 
reasons, the dispersion at the BPM location has to stay unchanged. But, when we 
optimize a cell in which periodic dispersion at the BPM location is by itself function 
of the cell phase advance, the situation, of course, changes. Let us, like in [1], con- 
sider a thin lens FODO cell of the length L in which two identical thin lens dipoles 
with transfer matrix 
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Figure 3: Functions QnifJ'p) shown for n = 3,4, 5,6 (magenta, red, green and blue 
curves respectively). 




(113) 



are inserted in the middle of drift spaces separating the focusing and defocusing 
lenses. Let us assume that the BPM is placed in the "center" of the focusing lens 
with the periodic dispersion at this locations being 



Vx,p{si) = V+ = 

where ip is the cell deflection angle. 
In this situation we can write 



Lip 1 + 1 sin(yUp/2) 
1 sm\fip/2) 



(114) 



n \L(p J 



(115) 



where functions depend only on the cell phase advance /ip and are converging 
(from above) to the function 



sin^(/ip/2) 



(l + lsin(/ip/2))^ 



:il6) 
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Figure 4: Functions shown for n = 3,4, 5,6 (magenta, red, green and blue 

curves respectively). The gray curve shows function \E'oo(/^p)- 

as n goes to infinity. 

The functions ifJ'p) for n = 3, 4, 5, 6 are plotted in figure 4 together with their 
values in the points (11121) shown as small circles at the corresponding curves. One 
sees that, again like in [1], there is nothing really special about points (I112p except 
the trivial fact that all of them belong to the graph of the function "^^o ■ 

Before switching to the situation when we have two BPMs per cell let us rewrite 
expression (llOSp for the function QnifJ^p) in the form 



/3p(si 



-1 



where 



1 

2\ 2 



is the mismatch between the error and the periodic betatron functions (even so we 
do not assume, in general, periodic betatron functions and/or periodic dispersion 
being the design betatron functions and/or design dispersion matched to our beam 
line) and 
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Figure 5: Functions <l>„(/ip) shown for = 3,4,5,6 (magenta, red, green and blue 
curves respectively). The gray curve shows function $oo(/^p)- 



'^Vl,pisi)ml{l3„ I3p) fl sm{nfip/2) 



n sin(yUp/2) 



1 sin(?7,/ip 



n sin(/ip) 



(119) 



is the difference between periodic and error dispersions measured by using the Courant- 
Snyder invariant formed out of error Twiss parameters. 

Note, for completeness, that if one will express the difference between periodic 
and error dispersions using Courant- Snyder invariant formed using periodic Twiss 
parameters, then one will have the following relation 

Ix {Ppy Vx,i; Vx,py Vps Vp,p) 



nip (/?,, f3p 



l_ sin{nfip) \ V -V ) 

n sin(/i ) / '^'P' 'P'« iP:PJ 



1 sm{nfip) 
in(/ip) 



n sm 



^ ^ 1 smjnfip) 
\ n sin(/ip) 



■ Ix (/3. 



(120) 
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Let us now turn to the situation when we have two BPMs per cell with 6 being 
the phase shift between the first and second BPM location. 

In this situation the square of the error energy spread can be expressed as 



A? 



bpm 



(121) 



where multiplier 



■U7n = \ I 



has a form which is very similar to flll7p with 



h (Z??, - Vx,p, Vp,'; - Vp,p) (122) 



4 [vlpisi) + rjl,p{s2)) ml /3p) / 1 sin(r2/ip/2) 



Pp{si) + /3p(s2) 



n sin(/ip/2) 



^ _ ]_ sin(ra/ip) 



n sin(/ip) 



Pp{si)ri%{si) + 2cos(g)v//3p(gi)/3p(g2)r/x,p(gi)r/x,p(g2) + /^p(g2) r/^,p(52) 
(/3p(si) + /3p(s2)) (r/2^p(si) + Tllp{s2)) 



+ 



+ 2sin2(^) 



/3p(si)/3p(s2) (I sm{nfip) 



(/3p(si) + /3p(s2)) sin(/^. 



(123) 



and with the mismatch between the error and the periodic betatron functions having 
now the following form 



rup {f3,, (3p 



1-1 - 4sin2(^) 



(3p{si)(3p{s2) \ f l_ sm{nnp 
(/3p(si)+/3p(s2))V \n' sin(/ip) 



2\ 2 



(124) 



For a thin lens FODO cell with BPMs placed in the "centers" of focusing and 
defocusing lenses we have 6 = /ip/2 and the periodic beta function and the periodic 
dispersion at the BPM locations are equal to 



1 ± sin(/ip/2) Llp 1 ± isin(/ip/2) 

P± = J^- —f — ^ , V± = —r 



sin(yUp 



sin (Aip/2) 



(125) 
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Figure 6: Ratio \[^oo / '^'oo as a function of the phase advance fip. 



With these assumptions we can write 

^bpm ( 4 



with functions converging to the function 

sin^(/ip/: 
2 + isin2(/ip/2) 



^ / N sin^(ur,/2) 

•^-(/ip) = . , 1 -t; ^o^ (127) 



as n goes to infinity. 

The functions $„ (/ip) for n = 3,4, 5,6 are plotted in figure 5 and one can see 
that though we are using two times larger number of BPMs, the energy resolution 
improves mainly in the region of the low phase advances, while for the high phase 
advances it stays almost unchanged. To understand the situation better, it is useful 
to look at the figure 6 where the ratio of the limiting functions \l/oo and $oo is 
shown. 



6 Courant-Snyder Invariant as Error Estimator 

When we consider the jitter problem, the subject of our real interest is the actual 
difference 8zq between parameters of the instantaneous and the golden trajectory. 
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Our measurement system delivers us an estimate of this parameter, which in- 
cludes the effect of the BPM reading errors ?. 

Thus, in the framework of the model considered, the only information which we 
could obtain about the true difference 6zq is the statistical information based on the 
properties of the random variable 6z^ — Szq, which, due to our assumptions, has zero 
mean and whose statistical distribution does not depend on the actual value of 6zo . 

It seems to be natural to use the module \6e^ — fcol as a numerical measure 
of the difference between estimated and true beam energies, but the quantitative 
measure of the difference 6z^ — 6zo from zero in the transverse phase space could be 
chosen differently. One may simply use the Euclidean vector norm, but, as it was 
already stated in [1], the usage of the Courant-Snyder quadratic form has certain 
advantages. For example, when one considers errors only in the reconstruction of 
the transverse orbit parameters in the beam line without dispersion, the Courant- 
Snyder quadratic form is an invariant and therefore all estimates based on it do not 
depend on the position of the reconstruction point. And, as one will see below, even 
for dispersive particle motion the Courant-Snyder quadratic form is a "much better 
conserved quantity" than the Euclidean norm. 

6.1 Transverse Jitter 

Let us first return to the situation whose study was already started in paper [1], 
where we considered errors in the reconstruction of transverse orbit parameters in 
the beam line without dispersion. 

Let /3o{r), ao(r) and 70 (r) be the design betatron functions, and 

I^{r,x,p) = 7o(r)x^ + 2ao{r)xp + f3o{r) (128) 

be the corresponding Courant-Snyder quadratic form. 

According to the above discussion, the object of our current interest is the random 
variable 

= 4(r, Sx^ - Sxo, 5p^ - 5po). (129) 
The mean value of this random variable was already calculated in [1] and is equal 

(4^> = 2e,mp(/3„/3o), (130) 

where mp(/3<^, (3q) is the mismatch between the error and the design betatron func- 
tions. That is, probably, all what one can obtain without making additional assump- 
tions about distribution of BPM reading errors. 

In this subsection we will assume that the random vector ? has a multivariate 
normal distribution and will find not only higher order moments of the random 
variable , but also its probability density. 

Calculations made in [1] show that we can represent the variable in the form 
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(131) 



where 



K = tv,m^r;\ 



"o/a/A) a/A) 



M 



( 



\ 



(Xri. Or, 



(132) 



and the components of the vector r/ = R~~^q are independent standard normal 
variables. The matrix K^K is by n matrix, but, as it was also shown in [1], it has 
only two nonzero eigenvalues, namely 



rrip (A, A) ± Y/m2(/3„/3o) - l) . 



(133) 



If e± are the unit orthogonal eigenvectors of the symmetric matrix K cor- 
responding to its nonzero eigenvalues ^± , then we can rewrite (11311) in the form 



(134) 



where C± = ^ are two independent standard normal variables. 

With representation f ll34p calculation of all probabilistic characteristics of the 
random variable becomes rather straightforward. For example, the following 
formula gives its variance 



= { {ID" > - (11)" = 4e,^ (2m5(/3„ /3o) - l) 



(135) 



Moreover, it is not very complicated to calculate the probability density of this 
random variable using, for example, results published in [6]. This density p{t) is 
equal to zero for negative values of its argument, and for t > 



2e, 



pit) = 7;-Io{ \/"^»(A, A) - 1 7^ exp -rUpilS^, /3o) — 



2e. 



2e. 



(136) 



where Iq is the modified Bessel function of zero order. 

Note that for mp{[3^, (3o) = 1 the density (11361) becomes the density of chi-square 
distribution with two degrees of freedom and in this case the distribution function 
F{t) can be calculated in the explicit form 



F{t) = Pr(4^ < t] 



p{T)dT = 1 - exp ( 



(137) 
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6.2 Transverse and Energy Jitter 

When the beam energy is included in both, measurement and dynamics, the 
transverse motion could be separated into two parts: dispersive motion and pure 
betatron oscillations. One can write 

6X0 = {6X0 - fco ■ Vx,0) + ■ 'nx,0 /-j^ggN 



and 

( Sx^ = {5x, - 5e^ ■ Tj^fi) + 5e^ ■ rj^fi ^^^^^ 



5p^ = (Sp^ - Se^ ■ ?7p,o) + fe? ■ Vp,o 

where rj^^o and rjp Q are the coordinate and momentum design dispersions respec- 
tively. 

The first terms in the right hand sides of formulas f ll38p and f ll39p represent the 
pure betatron oscillations. Let us at the beginning estimate their difference using 
the Courant- Snyder quadratic form, which in this case will be an invariant, i.e. let 
us consider the random variable 

= 4( r, {5x^ - 5xo) - {5e^ - feo) ■ Vxfi , (Sp^ - Spo) - {Se^ - Ssq) ■ ?7p,o )• (140) 

The mean value of this variable is given below 

(I'x) = 2 mp(/3^, /3o) + ■ 4 ( r, rj^^^ - r]^^ , rjp^^ - r]p^o ) , (141) 

and one sees that, in addition to the mismatch between error and design betatron 
functions, the difference between error and design dispersions starts to play an im- 
portant role. 

If we again will assume that the random vector ? has a multivariate normal 
distribution, we can represent in the form 

II = + (142) 

which is similar to fll34p and in which ^± are again two independent standard nor- 
mal variables. Unfortunately, the expressions for li± become essentially more com- 
plicated than fll33p and, with the notations 

rhp = mp{l3^, /3o), 4 = 4 ( r, r]^,^ - r]^fi , Vps - Vp,o ) , (143) 
are given below 



/i± = e,m, + ^ ■ 4. ±^62 {ml - l) + (m, - 1) . 4 + ^ . /2 . (144) 
With representation (11421) one can calculate the variance 
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(145) 



and also find formula for the probability density p(t). This density is equal to zero 
for negative values of its argument, and for t > 




(146) 



where Jq is the modified Bessel function of zero order and 

A = B = ~^±±^ (147) 

are the geometric and the arithmetic means of the eigenvalues fll44p respectively. 

To finish this section, let us note that in order to get probabilistic characteristic 
of the random variable fll29p . i.e. in order to study not the difference in the pure 
betatron oscillations, but the total difference in the transverse motion, one simply 
has to set to zero the design dispersions in all formulas of this subsection (indepen- 
dently, if actual design coordinate and momentum dispersions are equal to zero or 
not). The obtained formulas will, of course, not have invariant character anymore. 
Nevertheless, the dependence form the position of the reconstruction point will enter 
them through the single parameter, namely through the value Ix ( r, rjx^c , ''7p,<; ). 
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